99%-fidelity ballistic quantum-state transfer through long uniform channels 
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Quantum-state transfer with fidelity higher than 0.99 can be achieved in the baUistic regime of an 
arbitrarily long one-dimensional chain with uniform nearest-neighbor interaction, except for the two 
pairs of mirror symmetric extremal bonds, say x (first and last) and y (second and last-but-one). 
These have to be roughly tuned to suitable values x ~ 2N~^'^ and y ~ 2^''^N~^'^ , N being the chain 
length. The general framework can describe the end-to-end response in different models, such as 
fermion or boson hopping models and XX spin chains. 
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I. INTRODUCTION 

Transferring quantum states between distant registers 
is one of the basic tasks that a quantum computer based 
on qubits located on fixed positions has to accomplish. 
The general scheme where a quantum channel physically 
connects the sending and the receiving qubits naturally 
emerges from such requirement and a variety of propos- 
als have been recently put forward for realizing quantum 
channels for state transmission by different physical so- 
lutions: phonon modes for trapped ions 1] or photon 
modes for superconducting qubits [2, 3], molecular lat- 
tices for vibrational excitons ^4] or optical lattices for cold 
atoms [5-^ , arrays of coupled quantum dots [9l-[l2| , and 
interacting S — ^ spins on a one-dimensional lattice, usu- 
ally referred to as spin chains. In particular, proposals 
based on spi n chains were first introduced in this context 
by Bose [ij, |lj| and have attracted much attention in 
the last decade, due to both the possibility of exploiting 
their natural dynamics for the transfer process, and the 
availability of analytical results that allow for a detailed 
description of their dynamics. 

All of the above mentioned proposals are based on the 
idea that the state of the sender-qubit be transferred at 
large distance via a dynamical site-to-site hopping mech- 
anism that we will here describe in terms of a general 
hopping model. Once a 2-dimensional Hilbert space, gen- 
erated by |0) and |1), is assigned to each site i — 1, ...,N 
of a one-dimensional lattice, the hopping model is defined 
by the following Hamiltonian: 



1 ^ 



(1) 



where \i) = |0)®*~i|l)|0)®^"' is a brief notation for 
single-excitation states, and the structure of the hopping- 
amplitude matrix A — {Aij} depends on the properties of 
the specific model. Despite being responsible of a possi- 
bly long-distance transfer process, such mechanism can 



be generated by a short-range interaction. Notice that 
Eq. ([T]) describes free bosons as well as free fermions and, 
as far as the latter are concerned, an exact mapping ex- 
ists between the Hamiltonian ([T|) and that of the S — ^ 
XY chain, where the matrix elements Aij correspond to 
the exchange couplings. By the hopping model the basic 
transmission mechanisms can be studied and character- 
ized: this, in turn, can also enlighten on the specific role 
of nonlinearity J15l - [l7t , noise ^Ig] and dissipation jlg, |20| . 

The quality of the channel, as far as the transfer pro- 
cess is concerned, essentially depends on the hopping am- 
plitudes Aij] in particular, it has been proven that locally 
engineered hopping amplitudes can lead to perfect state 
transfer |2ll - [27f . However, strategies based on a very 
detailed design of the internal couplings appear very de- 
manding from a practical point of view 28| . 

Having in mind that experimental setups generally re- 
quire the couplings to be as uniform as possible (see, 
e.g., Ref. [29]), different authors have proposed alterna- 
tive strategies for obtaining high-quality state-transfer 
processes through mirror-symmetric channels with uni- 
form bulk and just a few extremal couplings allowed to 
be varied i3Qn35j . Some of these alternative strategies 
are based on the idea of markedly weakening the coupling 
between the channel and the sender/receiver qubits, a so- 
lution that has been shown 30, ^3^, |3g to lead to a very 
high-quality state transfer. However, this scheme yields 
large transmission times and requires such a reduction of 
the extremal bonds that a severe limitation on the ac- 
tual length of a functioning channel must be taken into 
account. 

In general, the transmission quality of almost uniform 
channels is expected to deteriorate as the length of the 
channel is increased [13, [ij, [3J| due to dispersion, which 
is integral to uniformly distributed couplings. On the 
other hand, modeling a scalable quantum-state transfer 
process whose quality depends as little as possible on 
the physical length of the channel is an essential issue to 
address, especially if solid-state implementations and/or 



experimental analysis are in order. 

Therefore, we find ourselves squeezed between the 
seemingly incompatible requirements of avoiding too 
much a detailed design of the physical channel and yet 
getting a reasonably long channel characterized by a rel- 
atively convenient transfer time. 

Our approach for solving this puzzle stems from the 
idea of exploiting, in an almost uniform channel, the bal- 
listic state-transfer mechanism that allows perfect trans- 
mission of a wavepacket |33l . 1371 . |38|, in virtue of a 
perfectly coherent, non dispersive, dynamics. Perfect- 
transfer follows from the requirement that the normal 
modes correspond to equally spaced eigenfrequencies (or, 
loosely speaking, to a linear dispersion relation), so that 
a coherent mirrored reconstruction of all normal compo- 
nents occurs [39J . We infer that only the modes involved 
in the initial configuration of the overall system need sat- 
isfying the linearity condition [33] in order to get an effec- 
tively ballistic dynamics: systems realizing this condition 
can be dubbed 'optimal state-transfer' systems. 

This idea has been first proposed and implemented |33l . 
kOJ I in a scheme where we could only play with one param- 
eter of the Hamiltonian, namely the value of the extremal 
bond X = An — Ajv_i,Ar. 

We have learned that the emergence of an optimal 
value for x follows from quite a complex interplay be- 
tween two conflicting effects, namely the deformation of 
the eigenvalue spacings and the shrinking of the mode 
distribution, which are simultaneously driven by the 
value of X. In order to further improve our results, we 
understand that these effects must be handled indepen- 
dently, a goal that can be accomplished by introducing 
just one more parameter in the model, as shown in this 
paper. Notice that the introduction of a second param- 
eter in the model should not be thought of as a way of 
moving towards the perfect-transfer scheme (which re- 
quires N j2 parameters) but rather as a practical answer 
to the effective need of controlling two competing effects. 

In Section |TT] the state-transfer mechanism in the hop- 
ping model Eq. ([T]) is studied and the quantities that 
characterize the efficiency of the quantum channel in 
terms of the transition amplitudes are obtained. In Sec- 
tion |llT] we describe the setup for implementing the two- 
parameter optimal-transfer scheme (with two adjacent 
modified bonds at the ends of the chain) and derive im- 
plicit analytical expressions for the frequencies and the 
matrix elements entering the transition amplitude. The 
behavior of the latter is thoroughly discussed in Sec- 
tion IIVI where numerical results are reported together 
with the analytical derivation of the large- A^ limit of the 
attainable optimal transmission fidelity. In Section |V] 
the dynamical behavior is shown to be ballistic, i.e., the 
information is carried by a wavepacket traveling at con- 
stant speed along the channel. Conclusions are drawn in 
Section rvTl Some details of the calculations are reported 
in Appendix El 



II. STATE TRANSFER IN THE HOPPING 
MODEL 

Let us consider the Hamiltonian ([IJ. In the case of 
nearest-neighbor interaction the matrix A is tridiagonal 
and, for the chain to be connected, the off-diagonal ele- 
ments cannot vanish. Without loss of generality [41] we 
assume A to be real with Ai^i^i — Ai^i^i > 0. Introduc- 
ing the orthogonal matrix U = {Uni} that diagonalizes 
the matrix A, 



ij 

and the one-excitation states 

i 

the Hamiltonian turns into the diagonal form 

N 

n = ^UJn\n){n\ . 

ri = l 



(2) 



(3) 



(4) 



If A is mirror-symmetric, i.e., it commutes with the 
mirroring matrix J = \Jij—^N+i-i,j\i then also J is 
diagonalized by [/, 



/ , ■^ij ^nj 



Un,N 



+ 1-1 



Jn ^n 



(5) 



and since J^ = 1 the eigenvalues j„ are either 1 or —1. It 
is proven in Ref. [42 | that if the eigenvalues LUn are cast 
in decreasing order then j„ = (— )"^^ = — e""; in other 
words, the n eigenvector of A is mirror-symmetric or 
-antisymmetric according to whether n is odd or even. 



A. Estimating state-transfer quality 

The purpose of state transfer is to start with a product 
state li/)) — \i{j)\0)'^^~^, with a generic state 



|V) = a|0)-f/?|l) 



(6) 



of the first qubit, and let it evolve with H in such a way 
that at a given time i ~ A^ the state of the last qubit is 
as close as possible to \ip). Formally, the evolved overall 
state is p{t) — e~*^*|V')(V'|e*^* and the instantaneous 
state at site i is 

Pi{t) = Tr{i_..._,_i_j+i^...^Ar}p(i) . (7) 

Defining |0) = 10)*^^, the evolved overall state obeys 

-»«t|„/,\ _ .Jn\ , «\ ^/"i^-^-w'IDIz) , (8) 



= a|0)+/3 V(*|e-*^*i 



so that 



P^{t) = 









(9) 



where Ui (t) is the transition amplitude from site 1 to site 
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Ui(t) 
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The state of the last qubit is described by the density 
matrix Pf^{t), whose degree of similarity with the initial 
state IV') (V' I can be estimated by the fidelity, 



F{^P,t) = {^|;\p,{tm 



(11) 



The overall channel quality can be estimated by taking 
the average of F{il;,t) over all possible initial states \ip), 
i.e., over the Bloch sphere, |ap + |/3p — 1. This yields 
the average fidelity J"(i) in terms of Uj^(i), with a de- 
pendence on the phase of the amplitude that one can get 
rid of by prescribing a proper local rotation to the last 
qubit [40|, or, in a s pin- chain context, by applying a suit- 
able magnetic field [lj|, in such a way that u„(t) is real 
and positive at the arrival time; with this proviso, setting 



L{t): 



At)\ , 



one has 



Tit) 



1 , [l+ujt)? 
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(12) 



(13) 



Hence, the average fidelity is a monotonic function of 
(the modulus of) the transition amplitude u{t). On the 
same footing, it can be shown that also the entanglement 
fidelity, that measures the efficiency of transfer for an 
entangled state with an external noninteracting qubit, 
reads lii 



-FE(t) = 



[l+u{t)] 



(14) 



It is clear at this point that u{t) E [0, 1] can be assumed 
as a transfer-quality indicator, u{t) = 1 corresponding to 
perfect transfer. For almost perfect transfer, u = l — e, 
one has a reduction of the fidelity to the first order in e, 



J^ ~ 1 — ■^ £ , J^E ~ 1 



(15) 



In the following we will thus concentrate onto the task of 
maximizing the transition amplitude u{t) for the partic- 
ular model introduced in the next section. 



B. Mirror symmetric case 

If A is mirror-symmetric, Eq. ([5]) tells that UnN = — 
e'^^'^C/ni, so the transition amplitude we are interested in 
takes the form 
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FIG. 1. Mirror-symmetric quantum channel consisting of A'' 
sites with uniform bonds of unit strength, but for two pairs 
of weaker bonds at the extrema, x and y. State-transmission 
is meant to occur between sites 1 and A^. 



which consists of a sum of time-evolving phases with 
a weight given by the mode distribution (or density) 
Vn- In the case of equally spaced eigenfrequencies a;„ = 
uJo + v{7r/N)n, where w is a constant, one would have 
u{t—N/v) = 1, i.e., full transmission of the state. As al- 
ready mentioned in the Introduction, this case of a 'linear 
dispersion' can be realized in discrete finite systems [2j] 
by allowing any single nearest-neighbor coupling to be 
properly tuned such that ^i,i+i oc ^Ji[N~i). 



III. THE QUASI-UNIFORM CHANNEL 

We assume from now on that the mirror-symmetric 
matrix A = {Aij} describes the chain of Fig. [U which 
is uniform except for the two extremal bonds A12 — 
An-i,n = X and A23 = ^Ar-2,Af-i = y, while all other 
couplings are equal and set the energy scale, Ai,i+i = 1 
for i = 3, ..., N—S. Such quasi-uniform interactions were 
used in Refs. [3, E^l for the purpose of studying long- 
distance entanglement in the ground state. Explicitly, 
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It is assumed that 0<a;<l and 0<y<l. In Ap- 
pendix]^ the A'' energy eigenvalues are found in the form 
LOn = COS kn , where the allowed values k = kn are deter- 
mined by Eqs. (|A12[) and (|A13p . and the mode distribu- 
tion Vk is given by Eq. (IA23[) . As the latter is symmetric 
around 7r/2, it is more convenient to work with the shifted 
variable 
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(18) 



(19) 



the corresponding eigenfrequen- 



cies are ujq ~ sing. Understanding that the variable q as- 



sumes the allowed values (H 
of the eigenvalue index n 



, we can use it in the place 
A^ -1-1/2 - m. Note that for 




0.4 , 0.5 
q/7T 



FIG. 2. (Color online) Shifts (ys,, Eq. ^U^, for different values 
of X and y. The dashed curves correspond to y = l. The fixed 
point gp — sin~^^ lies on the dashed straight line. 



even N the m's are half-integer: there is no qualitative 
difference between the outcomes for even and odd N, the 
latter case is just more easily handled numerically. The 
phase shifts 



"Pq 



tan 



y"^ sin 2(7 



2;2_(2-j/2)(i_cos2(7) 



2q (20) 



displace the g-values from the equally spaced values of 
the fully uniform case {(fq—Q for x = y = l). Note that 
the phase shifts do not alter the sequence of the g^'s, 
which increase with m. 

Note that the denominator of Eq. (pn|) can vanish and 
the argument of tan~^ can jump from oo to — oo, when 
this happens the conventional range of tan~^ has to be 
extended above 7r/2 rather than jumping from 7r/2 to 
—it/2. In this way Lpq is a continuous function, shown in 
Fig. [5] for positive q. The effect of weakening the bond 
y is evident: the deformation is such that Lpq change its 
convexity in an interval \q\ ^q-p{x) whose width is char- 
acterized by the fixed point q^{x) = sin~ |-. Therefore, 
the parameter y can change the separation between the 
allowed g-values, and, in turn, between the eigenfrequen- 
cies affecting the coherence in the time evolution. 

The mode distribution Vk^Ul^, Eq. (|A23p . can be 
written in terms of q as 

2 



N+l+2ip' 



(21) 
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X y 



v 

x^ + (4-a;2-22/2)2 tan^g - 16 (l-y^) sin^g ' 

The term in Lp' is clearly irrelevant for a long chain, even 
though it ensures the normalization, ^^^ Tq^ = 1, for any 
finite N . Fig. [3] describes the typical behavior of the 
mode density Vq when y is varied keeping a fixed x <1. 
Compared to its counterpart in the case of the uniform 
chain, x = y = 1, it is evidently more structured; by low- 
ering y its tails get increasingly suppressed and Vq def- 
initely changes to a bimodal shape with two symmetric 
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FIG. 3. (Color online) The normalized large- A'' mode density 
Vq, Eq. (|21[) . for a; = 0.6 and different values of the second 
bond coupling y. The thicker curves are those for y = l and 
for the threshold value y = Y(x) (see text and Fig. |4|. The 
broad dash-dotted line is the result for the fully uniform chain, 
x^y = l. 
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FIG. 4. (Color onhne) The function Y{x), Eq. ([22]), separates 
the occurrence of unimodal and bimodal mode distributions. 
The optimized pairs {x,y), i.e., those which maximize the 
arrival amplitude u{t), are also reported for selected values of 
N (see Table H). 



maxima: this occurs when in the denominator the coef- 
ficient of sin q becomes larger than that of tan^g, i.e., 
when y is smaller than 



Y{x)^dV2 



Y 



(22) 



which is the curve shown in Fig. 21 Vq is unimodal for 
y > Y{x) and has its maximum at g = 0; at the threshold 
y = Y(x) the maximum fiattens, the deviation being ~ 
g**, before developing the lateral maxima which are the 
more pronounced the smaller y. In the limit y— >0 the 
distribution Vq tends to two symmetric 5-peaks at g = 
±sin~^|-, corresponding to the excitations of the single 
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FIG. 5. (Color online) End-to-end amplitude at arrival u(x^ y) 
as a function of x for selected values of y and TV = 501. The 
thicker curves correspond to y = 1 and to the optimized value 
2/ = 0.5439. 




FIG. 6. (Color online) Contour plot of the average fidelity at 
arrival F{x,y) in the {x,y) plane for A'^ = 51. 



dimer with interaction x, which is indeed isolated from 
the chain when y = 0. 

To represent the behavior of the level spacings, whose 
uniformity is crticial for the coherence of transmission, 
one can define a sort of 'group velocity' by 



N+1 



J drnl = 



N+\ 



7V-hl-h2(^: 



7 cosg , 



(23) 



where the 'derivative' of Eq. dH]), {'N^'\+2l^'^^m<l = 'l^■, 
has been used, and from Eq. (|20l) 



V', - -2 - 



2y2[a;2 + 2(2-a;2-2/2)sin^5] 



a;4 -F 4[y4_a;2(2-y2)] sin^g -^ 16(l-y2) sin^g ' 

(24) 




IV. TRANSITION AMPLITUDE 

A. Numerical results 

The transition amplitude between the sites 1 and iV at 
the time t is given by Eq. (IT51) . It can conveniently be 
rewritten as a sum over the allowed (/'s (J19p . 



^w = |E^^ 



^i{7n7i—t sin Qtii) 



(25) 



This sum can be evaluated numerically: for any pair 
(x, y) our code solves iteratively the coupled Eqs. P^ 
and (|20|) and looks for the value oi t>N, i.e., the ar- 
rival time, where u(t) attains its largest value, say u{x, y). 
Typical results for u(x, y) are reported as a function of x 
for selected values of y in Fig.jSJ which refers to a chain of 
length iV = 501. One can see that taking y smaller than 1 
the amplitude can become much closer to 1, also for the 
longest channels, and it is possible to identify the optimal 
values (x'^P*, y'^P*) that make u{x,y) reach its maximum 



FIG. 7. (Color online) Contour plot of the average fidelity at 
arrival J-{x, y) in the {x, y) plane for A'^ = 501. 
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Moreover, the maxima are so broad 



in the (x, y) plane that a relatively large mismatch from 
the optimal values still keeps giving a large amplitude, 
and consequently a large average transmission fidelity, as 
can be appreciated in the contour plots of Figs. [Hand [T] 
The numerically evaluated optimal data are reported 
in Table HI together with those obtained in Ref. [4^ by 
varying only x with y = 1. Comparing with the latter, it 
is seen that the transfer quality improves in an extraor- 
dinary way: even for a chain of length iV = 50 001 the 
amplitude increases from 0.859 to 0.987, i.e., the state- 
transfer fidelity F rises from 0.909 to 0.991. The values 
of Table |T] show that the optimized transmission ampli- 
tude decreases more and more weakly for large N , so that 
there could be a finite asymptotic amplitude attainable 
for an infinite chain: this is the case, indeed, as shown 
in the next subsection. In Fig. |5] the optimal values for 



TABLE I. Optimal pairs (2;°^', j/°''*) and the corresponding amplitude u°p' : 
channel lengths TV. Also reported are the optimal values obtained [40|]. 
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FIG. 8. (Color online) Optimal values for x, from Table [Jl 
reported vs. A'^ both for j/ = 1 and for the optimized y = y°'^^ . 
Note that a;~ A~^'® in the former case, while x^N~^'^ in 
the latter. 



X are reported for y = 1 and for y — j/°p* . As shown for- 
mally in the next subsection, it turns out that 2:°p* scales 
as N~^/^ in the former case, while it obeys a new scaling 
law, apparently N~^/^, in the latter. 

The great improvement in transmission quality deals 
with the same argument we gave in Ref . 4(j] : although a 
constant group velocity Vq yields perfect transmission, it 
is sufficient that Vq be constant for the (/-modes excited 
by the initialization of the first qubit whose distribution 
is Vq ■ The results illustrated above confirm that the pos- 
sibility of controlling, by means of two parameters {x and 
y), the two most relevant features, shape oiVq and stabil- 
ity of Vq , allows us to obtain an optimal trade-off leading 
to nearly perfect transmission. 

The weaker second bond y acts indeed on the group 
velocity, as it appears in Fig. [SI where the shape of Vq 
is reported for the optimized value of x, and compared 
with the corresponding mode density Vq. With a smaller 
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FIG. 9. (Color online) Group velocity Vq and corresponding 
mode density Vq for A = 501, x = x' 
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■ 0.225, and selected 



values of y. The thicker curves correspond to y = 1 (dotted) 
and to the optimized value y — j/"^* — 0.544 (solid). 



y, the central dip appearing for ?/ = 1 is strongly reduced 
and the group velocity just shows a small modulation in 
a rather wide range, so favoring a coherent dynamics; at 
the same time Vq broadens, but not dramatically as it 
is mainly controlled by x ( see, e.g.. Fig. 4 of Ref. [40{): 
the optimized value y = y°P' clearly gives the best result 
compatible with the assumed parametrization. 



B. Asymptotic behavior for large A'^ 

The numerically estimated optimal pairs (x°P',y°P') 
are shown in Fig. 2] They evidently lie almost exactly on 
the threshold curve Y{x) separating the unimodal from 
the doubly-peaked shape Pq, i.e., y°P* ~ y(a;°P*). We ex- 
pect this to hold in the large- TV regime, since the condi- 
tion y = Y(x) clears the quadratic terms in the denomi- 
nator of the density (PT|) leaving smaller tails ~ q~'^, thus 
allowing for a more effective cut of the modes involving 
the main nonlinear part of Vq. Therefore, instead of con- 
sidering y as a free parameter, we fix it to be given by 



2/ = r (:e) = W V2 a:-y = 2^/4 a;i/2 + 0(^3/2) . (26) 
As iV — > oo the sum ((25l) can be written as an integral, 

„i(irT?i-tsing„) /27) 



Uooit) — lim / dm Vn e" 

Af-i-oo ' ^ 



In order to evaluate this asymptotic expression, first note 
that from P^ one has to set 11171 = {N+l)qm -\- 2(/3g^, 
and TT dm ~ (iV+l+2i^' ) dq, so that 



Uoo{t) = lim 

Af^oc 



^ -p ^t[(N+l)q+2v^-tsinq] ^ /28) 
TT 



where 



Va^ 



2a;2y2 



x^ + {A—x'^—2y'^y tan^q — 16 (1— y^) sin q 



(29) 



is exactly the normalized function reported in Fig. [3l As 
by increasing N the optimal distribution gets narrower 
and narrower, the denominator can be expanded taking 
into account the assumption (j26|) : 



Va 



23/2^3 
x4 + {2qY 



(30) 



hence, the width of the relevant g-region shrinks with x. 
Let us introduce the scaled variable S, = 'i.q/x, which is of 
the order of unity, so that 



Vq dq' 



V2d^ 

1+e 



(31) 



As for the phase in Eq. ([28| . the leading term of the 
expansion of Eq. (^01) is 



(fiq ~ tan 



1-e 



(32) 



and, defining the arrival-time delay s by t = N+l+s, the 
remaining terms read 

{N+l)q - ts'mq = t{q - s'mq) -~sq~T^^~a^, (33) 

where 

1 /a;\3 



6 



(I)' 



t 



X 

a = — s , 



(34) 



are the rescaled counterparts of the arrival time t^ N 
and delay s^N^^^. Eventually, the asymptotic value of 
the amplitude reads 

V2 7,^ exp[z(T^3_^^+2tan-i^)] 



Woo(t,ct) = dS, 



1 + ^4 



(35) 
For a numerical evaluation it is convenient to perform the 
substitution ^ = tanz, and consider the maximization of 



2V^ p l-Htan^z 
Uoo = / dz ^ ^ ^ ^4 cos$(z) , (36) 



TT 



l+tan''z 



with 



<i>(z; T, cr) = r tan z — a tan z + 2 tan 



tan2z 

~72~ 



; (37) 



although this phase strongly oscillates for z close 
to 7r/2, the weighting function makes the numeri- 
cal convergence easy. The overall maximum corre- 
sponds to (T,cr) = (0.15545, 3.1645), and amounts to 
Woo — 0.987153, which is the asymptotic value reported 
in Table H] together with the asymptotic scaling resulting 
fromEq. (jM)) . 



/fi-r\ 1/3 

^opt^2(^) ~ 1.954 iV-i/3 



y 



'°pt = r(a;°P') ~ 1.662 N^^^^ , 



(38) 



while the delay scales as s = 2o'/a;~ 3.239 N^^^, so that 
the arrival time is 



t~7V + 1 + 3.239 N^/^ 



(39) 



STATE-TRANSFER DYNAMICS 



In Section |lT] we introduced the instantaneous transi- 
tion amplitude from site 1 to any site i of the chain, 
Eq. pO|) . This quantity substantially tells where the in- 
formation concerning the initial quantum state sits at 
any time t. Indeed, it obeys the sum rule 



N 



J2\uM'-Y.U'nl-l 



(40) 



n=l 



SO one can view the state-transfer process as the transmis- 
sion of a traveling wavepacket of amplitude Ui (t) which is 
able to optimally rebuild most of its content in the single 
A^-th site. Looking at the dynamics of this wavepacket 
sheds further light onto the ballistic transfer mechanism 
and allows for interesting comparisons, based on the data 
for \ui{t)\ reported in Figs. [TOl and [TT] 

The first comparison is made in panels (a) , (b) and (c) 
of Figs. [TO] and [Til and involves: (a) the fully uniform 
channel, which displays a very dispersive dynamics and 
is indeed inefficient for transmission; (b) the channel 



with only one optimized extremal bond [40J, showing 
an increased coherence; (c) the further step with two 
optimized extremal couplings which improves transmis- 
sion close to perfection. These features can be appre- 
ciated looking at the height of the arrival-time maxima 
in Fig. [ini and at the the rugged features, more evident 
in Fig. [TI] that represent the amplitude losses due to 
dispersion, according to Eq. (HUI . Fig. [11] makes also ev- 
ident the increasing arrival delay from (a) to (b) and to 
(c), as a consequence of the slower packet injection and 
reconstruction due to the softened endpoint couplings. 

Eventually, let us consider the case of perfect state 
transfer [2j], which is obtained by designing all nearest- 
neighbor couplings along the chain proportionally to the 
height of a semicircle of radius N+1 drawn over the 
chain. 



A, 



i,i+l 



= A 



i+lA 



N- 



-VW^) 



(41) 



the energy unit being arbitrary, it is chosen here such 
that the resulting linear spectrum 



iV-1 



N-1 



iV-l-1 



(42) 



yields the exact arrival time t — N+1. Notice that, ne- 
glecting terms ^ N^^, the maximum of the couplings 
{Ji} is 7r/2, their average is 7r^/8, and their mean-square 
value is 7r/3; this allows for a meaningful comparison with 
our model, and the corresponding data are reported in 
the panels (d) of Figs. [TU] and [TTJ The modulated cou- 
plings determine a varying velocity of the wavepacket 
along the chain, at variance with its constant velocity 
in the uniform channel of panels (a-c) . 



VI. CONCLUSIONS 

We have shown that almost perfect ballistic quantum- 
state transfer with fidelity larger than 0.99 can be ob- 
tained in an unmodulated channel of arbitrary length 
TV, just by allowing the two endpoint pairs of nearest- 
neighbor interactions, x and y, to assume optimal val- 
ues (a;°P*, 2/°P*). In addition, this maximum of the trans- 
mission quality in the (x, y) plane, as measured by the 
average fidelity or equivalently by the transition ampli- 
tude (jl6|) . is so broad that an experimental realization 
would not be bound to a fine-tuning of the endpoint cou- 
plings, as Figs, inland [7] clearly show. 

One might think that the approach presented here 
looks like the second step of a sequence that, by allow- 
ing further bonds to vary, would lead to perfect ballistic 
transmission as found in Ref. [23|. However, such de- 
duction does not hold true: indeed, while the aim of 
Ref. [23] is that of obtaining perfect transfer by letting 
all normal modes to evolve coherently, here we look for 
high-quality transfer by requiring that only the modes 
excited by the initialization of the first qubit be able to 
evolve coherently. This is done within an effective scheme 




FIG. 10. (Color online) Space-time perspective views of the 
propagating wavepacket |ui(f)|, Eq. (|10p . for a chain of length 
N = 251, (a) in the fully uniform chain x = y = l, (b) in the 
case J/ = 1 and optimal x = 0.422 [4Q|, (c) in the quasi-uniform 
channel, Eq. UTT} . with the optimal a:: = 0.276 and j/ = 0.598, 
and (d) in the perfect-transfer channel |23|], Eq. (|4ip . 
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ruled by two effects, namely the mode distribution and 
the frequency spacings, that can be kept under control by 
just two parameters, i.e., the extremal couplings x and y. 
Our approach leads to extremal couplings a;~ N~^^^ and 
yr^ N~^'^, sensibly larger than those required for perfect 
transfer, Eq. (|4T1) . which scale as N^^/"^. 

The asymptotic results improve from the uniform chain 
{x^y^l) which gives fidelity J^oo = 1/2, to the single 
tuned bond [40| (y = l) giving J^oo = 0.902, and eventu- 
ally to the optimal asymptotic fidelity evaluated in this 
paper, J"oo = 0.991. Whether comparable results could 
be obtained by other setups is difficult to say: the vari- 
ants are numerous and we cannot tell a conclusive word 
in this respect. 

The quasi-uniform channel here considered was previ- 
ously used in Refs. @, |43| for the different purpose of 
exploiting the quasi-long-distance entanglement shared 
by the extremal spins in the ground state in order to get 
efficient teleportation; in those papers, proposals for re- 
alizing the model by means of coupled cavity arrays and 
ultracold atoms in ID optical lattices are put forward. 
We expect that such physical realizations might possibly 
fit also our setup. To this respect, we remark that when 
the presence of disorder in the couplings along the chain 
is taken into account, the one-parameter optimal state- 
transfer scheme (y = 1) results to be more robust than 
the perfect-transfer one |J_8|; we expect such advantage 
to be preserved in the two-parameter optimal-transfer 
scheme. 

We finally mention the similarity between the problem 
of quantum-state transfer and the subject of continuous- 
time quantum walks [4J| , where regular space-time struc- 
tures called quantum carpets |45| can emerge from a 
complex dynamics: for instance, the revival of the wave- 
function is analogous, in a mirror-symmetric context, to 
state transmission. 



FIG. 11. (Color online) Space-time contour views of the 
propagating wavepacket |«i(i)|, Eq. (fTO)l in the same cases 
of Fig. [To] 
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Appendix A: Spectral solution 
1. Characteristic polynomial 

The characteristic polynomial x„(A;x,2/) = det[A 
A{x, y)] associated to the matrix pT)) is 

A ~x 
-X A -y 



X«(A) = 



-y A -1 
-1 A -1 



-1 A -y 

-y A ~x 
—x A 



(Al) 



N 



Expanding x^ (A; x, y) in the last column one finds 

X„ = (A^ - x^) C_, - Aj/2 C„^3 : (A2) 

where ^„ is the characteristic polynomial of the associ- 
ated matrix with only one nonuniform endpoint, 



C(A) 



A 


— X 












-X 


A 

-y 


~y 

A 

-1 


-1 
A 


-1 
-1 


A 

-1 


-1 
A 



(A3) 



N 

which can be expanded in the very same manner getting 
the analog of Eq. (fX2|) . 

C« = (A^ - x^) Vm-^ - Ay^ ?y„_3 , (A4) 

in terms of the characteristic polynomial rj^ (A) = 
X„(A; 1, 1) of the fully uniform matrix A^ (1, 1). For ?7„ 
one has the recursion relation 



^« = A 7?^ 



V„ 



(A5) 



together with the conditions -q^ =1 and 77^ = A, it can be 
solved in terms of Chebyshev polynomials of the second 
kind, 



^A 



sin(7V+l)A: 
sin k 



A = 2 cos k 



(A6) 



Hence, in the uniform case the N solutions of the secular 
equation Xiv(A; 1, 1) = VnW — correspond 46] to the 
following discrete values of k: 



TTTi 

N+l ' 



1, 



,N . 



To extract the dependence on x and y in Eq. 
first uses Eq. (j^H, 

X„(A) = {\^-x^)[{\^^x^)r,,_, - \y^ri,_ 
-Ay2[(A2-a;2)r;„_,-AyV- 



(A7) 
2]) one 

,(A8) 



and then Eq. (jA6p in the form rj^ sinfc — Sjje*'^'*'^-''^}, 
X„(A) sinfc = 3{(A2-.T2-Ay2e-''=)'e^(^-3)*} 



= 9f/„? AN+l)k 



K 



} 



(A9) 



where 



2. Shifts of the eigenvalues 

Eq. (jA9l) leads to a compact expression of the secu- 
lar equation, Gje**^^"''"'^^'^ u^.} = 0; in the fully uniform 
case, when both x,y — ^ 1, it has the solutions given by 
Eq. ISI- Setting 



Uk=\uk\e-'^'^ , 



(All) 



the secular equation reads 3{e*[(^+i)'=-2¥"=l} =0, and 
the eigenvalues A„ — 2 cos kn can be expressed as devia- 
tions from the uniform-case values (jA7|) due to the phase 
shifts ipk , 



kn — 



N+l 



n — I, 



,N . 



(A12) 



Keeping in mind that the variable k actually assumes 
the N discrete values {fc„}, one can unambiguously use 
the index k in the place of n. An explicit expression of 
the phase shifts follows immediately from Eq. (jAlOp . by 
separating the real and imaginary parts of Ufce^**^, 



ipk — 2k — tan 



y^ sin2fc 



(2-a;2-y2) + (2-y2)cos2fc 



Mode distribution 



(A13) 



To evaluate the transition amplitude (1161) one needs the 
square components of the first column of the orthogonal 
matrix U defined by Eq. ([2]). A nice formula derived 
in Ref. [4l[ (Corollary 7.9.1) comes to help: with the 
formalism used here it reads 



In — U^i 



(A14) 



9ax„(A„) 

where X2n(A) is the first minor of the determinant (jAl|) . 
It is equivalent to use the variable k= cos^^ ^ e (0,7r) 
and write 



Vk = C/|i = -2sin/c 



X2:iv(fc) 

dkXNik) 



(A15) 



with k definitely taking the allowed values (jA12| . i.e., 
Xjv(fc) = 0. As X2.iv(A) = A^„_2-y^Cjv-3 7 with calcula- 
tions similar to those which lead to Eq. (|A9p it is found 
that 



X,.^Ak)smk = ^{e''"'ukVk} 



(A16) 
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with 

Vk = l + (l-y')e-'''= . (A17) 

By deriving Eq. (|A9p with respect to k one has 

sinfc duxAk) = 3{e'(^+i)'= [i{N +l)ul + 2uku'^]) 
= (iV+l)5R{ 6^(^+1)'= M^} + 

+ 2 5{e*(^+i)'=7.fc4} ; (A18) 
the argument of 5R is real by the secular equation so, using 



sinfc 9fex„(fc)= [iV+l + 2 3f^| 
L Lwfc J 



,»(W+l)/c ^2 



(7V+l-2^1.)e'^^+'^'=ul. (A19) 



Eq. (JA15I) becomes 



^fc 



2 sinfc '^{e'^'^ukvu} 



2 sinfc '^{e^'^Ukvl} 



N^l-2^' 



\Uk\ 



(A20) 



By means of Eqs. (JAIOI) and (IA17P one can express 
Uk = Wfc(l+e~^*''') — x^e~^*'^, and a simple expression of 
the numerator follows, '^{e^^Ukvl} = x^y^ sinfc, finally 
yielding 



Vk = 



2a;^y2 



sin k 



N+1^2^'^ |ufc|2 ■ 



(A21) 



J 



A manageable expression for |ufep arises by working out 
Eq. dXlOl), 



e^-'^uk 



-x2+2(2-?/) cos^fc + 2i y^ sin fc cos fc , (A22) 



and taking the square modulus, after some algebra the 
outcome is 



Vk = 



2a;2y2 



iv+1-2^;^ 



(A23) 



x^ + (4-a;2-2y2)2 cot^fc - 16 (l-y^) cos^fc ' 

An explicit expression for ip'^ can be found going back 
to Eq. (lil9)) 



Vk 



Ui. 



»{<■<■} 



(A24) 



,Uk) \Uk\ 

where, from Eq. (IA10|) . 

4 = -2z[«fe-l-(lV)e-*'1 , (A25) 

and some further calculation gives 

Kl Vfe - 25RK K - 1 - (1 V)e-^^1 } 

= 2|Mfc|2 - 22/2[a;2 + 2(2-a;2-y2) cos^A;] . (A26) 

Eventually, the mode density can be made fully explicit 
starting again from Eq. (IA21[) . 



Vk = 



2x^y^ sin^fc 



{N-3){[x^ - 2(2-y2) cos2fc]2 + 42/4cos2fcsin2fc} + 4y2[a;2 + 2(2-j;2-y2) cos2fc] 

I 



(A27) 
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